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Abstract. For a Brownian system the evolution of thermodynamic quantities is a 
stochastic process. In particular, the work performed on a driven colloidal particle 
held in an optical trap changes for each realization of the experimental manipulation, 
even though the manipulation protocol remains unchanged. Nevertheless, the work 
distribution is governed by established laws. Here, we show how the measurement of 
the work distribution is influenced by the presence of rotational, i.e. nonconservative, 
radiation forces. Experiments on particles of different materials show that the 
rotational radiation forces, and therefore their effect on the work distributions, increase 
with the particle refractive index. 
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1. Introduction 

Most physical, chemical and biological phenomena present intrinsic degrees of 
randomness. Such randomness can be treated within the framework of statistical physics 
and thermodynamics. Equilibrium thermodynamics, in particular, is a mature discipline 
with well-established laws [HE]. However, most natural and engineered processes occur 
far from equilibrium, e.g., biomolecular reactions in molecular machines [3], oscillations 
of pumped mesoscopic chemical reaction systems [I] and fast switching between phonon 
distributions in optical cavities [5]. Such processes cannot be treated within the 
framework of classical equilibrium thermodynamics. This fact has motivated an intense 
activity during the last years aiming at extending the results of thermodynamics to 
out-of-equilibrium systems [61 [Tl El El HQl HU [121 H31 US UHl [16]. 

In this context, the study of thermodynamic quantities in microscopic systems, 
such as biomolecules and nanomachines, has posed new challenges due to the ineluctable 
presence of a Brownian noise background, which prevents one from straightforwardly 
scaling down the approaches employed on macroscopic systems. For example, the 
presence of Brownian noise alters the measurement of forces acting on microscopic 
objects, leading to artifacts if not correctly taken into account (T7J. Also the evolution 
of thermodynamics quantities, e.g., work and entropy, although deterministic in 
macroscopic systems, becomes stochastic in the microscopic realm. Indeed, a different 
value is measured in each experimental realization; even phenomena that are forbidden 
in a macroscopic systems become possible at a microscopic level, e.g., entropy-decreasing 
trajectories [T81 IT9"]. Nevertheless, the probability distributions of these thermodynamics 
quantities obey some deterministic laws [191 1201 IZD [221 [23] . 

Optically trapped particles [21] have emerged as a powerful model system to address 
experimentally novel concepts in the context of statistical physics in a convenient way 
relying both on the presence of a natural noisy background and on a finely controllable 
deterministic optical force field [25]. Indeed, a colloidal particle, i.e. a microscopic 
particle suspended in a fluid, is set in constant movement by the presence of thermal 
fluctuations, which are responsible for its diffusion and introduce a well-defined noisy 
background. Furthermore, it is possible to make use of optical forces to introduce 
deterministic perturbations acting on the particle in a controllable way. For example, 
optically trapped particles have been used to study experimentally many non-trivial, 
sometimes constructive, aspects of the presence of noise, e.g., stochastic resonance [26J, 
stochastic activation [27], Brownian ratchets [28], or stochastic resonant damping [29] . 

An implicit assumption in these studies is that the optical trap acts as a harmonic 
potential, i.e. the restoring optical force is proportional to the particle displacement 
from the trap center. This entails that the optical force field should be conservative, 
excluding the possibility of a rotational component. This is actually true to a great 
extent in the xy-pl&ne perpendicular to the beam propagation direction, at least for a 
standard optical trap generated by a Gaussian beam. However, due to the presence 
of the scattering force, this has been shown not to be true in a plane parallel to the 
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beam propagation, e.g., xz-plane (5UJ ED E2J EH]- For most experimental situations 
the nonconservative effects are effectively negligible, becoming significant only when the 
particle is allowed to explore a large region of the optical field [33] . This typically occurs 
for extremely low-power trapping, far away from the trapping regimes usually used in 
experiments, or for particles driven out-of-equilibrium. 

In this article, we show paradigmatically how the determination of a thermodynamic 
quantity, namely the work, is affected by the presence of a rotational, i.e. 
nonconservative, force field. We measure the work performed by a fluid flow on three 
optically trapped colloidal particles of comparable dimension, but different refractive 
index. For a larger refractive index, the rotational optical force field component increases 
and, therefore, introduces larger artifacts into the measured work distributions. 

2. Experimental setup 

The experimental setup consists of an optical tweezers build on a home-made optical 
microscope with a high-numerical-aperture water-immersion objective lens (Olympus, 
UPLAPO60XW3, NA=1.2); the optical trap is generated by a frequency and amplitude 
stabilized Nd-YAG laser (A = 1.064 /xm, 500 mW maximum output power, Innolight 
Mephisto)[M]. 

We use colloidal spheres of comparable diameter, i.e. ~ 1 /xm, of three different 
material - silica (Si), polystyrene (Ps) and melamine (Me). The main parameters 
of the experimental setup and of the beads used in this experiment are reported in 
table [TJ The particles are diluted in distilled deionized water to a final concentration 
of a few particles//il. The sample cell is made with a 150 /xm-thick coverslip and a 
microscope slide, which are separated by a 100 /im-thick Parafilm spacer and sealed with 
vacuum grease to prevent evaporation and contamination. Such sample cell is mounted 
on a closed-loop piezoelectric stage (Physik Instrumente PI-517.3CL), which allows 
movements with nanometer resolution. Obviously, a movement of the stage in a given 
direction corresponds to a movement the optical trap focus in the opposite direction 
from the perspective of the trapped particle. The sample temperature is continuously 
monitored using a calibrated NTC thermistor positioned on the top surface of the 
microscope slide and remained constant within 0.5 K during each set of measurements. 

A colloidal sphere is trapped and positioned in the middle of the sample cell, i.e. far 
away from the surfaces to avoid spurious effects arising from diffusion gradients [TT] 135] . 
Its x, y and z coordinates (figure GD^a)) are monitored through the forward scattered light 
imaged at the back focal plane of the condenser lens on a InGaAs Quadrant Photodiode 
(QPD, Hamamatsu G6849) [36j [37]. A digital oscilloscope (Tektronix TDS5034B) is 
used for data-acquisition. The QPD-response is linear for displacements up to 300 nm 
(2nm resolution, 250kHz bandwidth). The conversion factor from voltage to distance 
is calibrated using the power spectral density (PSD) method [35] . 

To minimize external noise, the experimental setup is mounted on a passive 
vibration isolation optical table. Moreover, the laser paths are kept as short as possible 
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Table 1. Particle and optical trap parameters. k p is the radial trap stiffness, k z is the 
radial trap stiffness, ry is the ratio k z /k p , rj their ratio and e the relative contribution 
of the nonconservative force fields [33J . 



Silica 



Polystyrene Melamine 



diameter (/im) 
refractive index 
density (g/cm 3 ) 
k p (pN//xm) 
k z (pN//zm) 



0.97±0.01 
1.37 
1.96 



1.38±0.09 
0.26±0.01 
0.19±0.01 



0.99±0.02 1.00±0.02 



1.17±0.07 0.84±0.08 

0.27±0.01 0.27±0.01 

0.23±0.01 0.32±0.01 

2.2±0.1 6.2±0.2 

0.4 ±0.02 0.4 ±0.02 



1.59 1.68 



1.06 1.51 



V 




0.5 ±0.02 



to avoid pointing fluctuations and the laser beam is enclosed in plastic pipes wherever 
possible. Finally, all the setup is enclosed by a polystyrene box to prevent air circulation 
and temperature drifts. In order to check the achieved stability, in figure [2] the PSD for 
the x (blue squares) and z (red diamonds) coordinates are presented for the Ps particle. 
They fit well the expected theoretical curves (lines) in the range between 50 mHz and 
500 Hz. In particular the PSDs show a very flat plateau up to the minimum frequency 
measured which corresponds to a period of 20 s. This allows us to exclude the presence 
of low frequency noise, since, as described above, each single measurement is performed 
within 10 s. 

3. Manipulation protocol 

The duration of each trajectory measurement is 10 s (figure HJc)). During the first 
5 seconds the stage is at rest and the particle is in an equilibrium state in the optical 
potential well, U(p) = k p p 2 /2 (figure [Qa)), where p is the radial position of the Brownian 
particle with respect to the center of the trap, i.e. p = yjx 2 ± y 2 , and k p is the radial 
stiffness of the optical trap, assumed to be harmonic. At t — 5 s the stage starts moving 
at a speed v = 1 pm/s along the x coordinate (figure [H(b)), producing an effective fluid 
flow and an effective force F = acting on the particle, where 7 is the particle friction 
coefficient. Hence, the manipulation protocol is 



After a pause of 1 s the above described sequence starts again, but the stage is moved 
in the opposite direction. This allows us to remain in the same region on the sample 
cell avoiding systematic drifts. By repeating this procedure 400 times, we acquire 800 
trajectories for each particle. 




t < 5s 
t > 5s 



(1) 
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During the application of the force, the average trajectory of the particle (black 
solid line in figure QJc)) climbs the optical potential; however, due to the overwhelming 
presence of the Brownian motion, each realization (e.g., gray solid line in figure [U^c)) 
presents a random behavior. 

Assuming a harmonic trapping potential, as the stage moves with constant velocity 
v along the x direction, the work done on the Brownian particle reads [19] 



W(t) = f dsvk x x{s), (2) 

Jt 

corresponding to the classical expression "displacement times force," where ds v is the 
infinitesimal displacement of the center of the potential and k x x(s) is the force acting 
on the particle with k x the trap stiffness along the x direction. 

In the case of an ideal Gaussian trap, the work probability distribution is given by 
HUH9] 

(3) 



$(W,t) = , 1 exp 



2a(ty 



which is a Gaussian whose mean and variance depend on time, i.e. /i(t) = 
v 2 r x k x (e~ t / Tx — 1 + t/r^j and cr(t) 2 = 2nkBT[i(t) with r x = j/k x the relaxation time in 
the optical trap, T the temperature and ks the Boltzmann constant. As will be shown, a 
deviation from a harmonic trap, leads to deviations from such work distribution. Thus, 
in the presence of nonconservative optical forces, we expect the experimental probability 
distribution function of the work to deviate from eq. [3j 



4. Work measurement on driven Brownian particles 

From the measured trajectories we calculated the work performed by the fluid flow on a 
Brownian particle according to eq. EJ The results are plotted as histograms in figure [3] 
for Si (first row), Ps (second row) and Me (third row) particles and for t = 10ms (first 
column), t = 500ms (second column) and t = 3000ms (third column). The theoretical 
expectation for a conservative force field according to eq. [2] are plotted as solid lines. For 
t — 10 ms, the work performed on all particles is in good agreement with the theoretical 
expectation. However, for larger times only the data relative to the Si particle keep on 
being in agreement, while deviations arise for the data on Ps and Me particles. 



5. Influence of rotational force fields 



In order to understand such disagreement, we consider the deviation of the optical force 
field from the harmonic one assumed in eq. El Indeed, in the case of optical trapping 
powers as small as the ones we employ, the optical force field can be better approximated 
allowing for a rotational component in the vertical plane |33j. 

A generic force field can be characterized by measuring the auto-correlation and 
cross-correlation functions between the coordinates of the particle [3E1 EHl HO]. In 
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particular, here in order to characterize the rotational optical forces, we use the cross 
correlation difference between p and z [33] 



7 l + k z /k p 



C pz (r) = 2 - ; - lu exp 



sinh 



f- 2 -€ 2 T7 2 \T 



P 



t — e t 



(4) 



where f = |[r p + r z ] is the average relaxation time in the trap, r z = j/k z the relaxation 
time along the ^-direction, k z the stiffness along the ^-direction and e the relative 
contribution of the nonconservative force fields. C pz {r) = is clear evidence that the 
contribution of the nonconservative force field component is effectively negligible. 

In figure 0] the experimental C pz (r) is presented for the case of a Si (blue), Ps 
(red) and Me (green) optically trapped Brownian particle without fluid flow. While the 
rotational component is negligible for the case of a Si particle, it becomes significant 
for the cases of Ps and, even more, for the case of a Me particle. Summarizing, the 
rotational component appears to increase with the refractive index of the particle. The 
same results are obtained in the presence of a fluid flow. 

The presence of a larger rotational component of the force field appears to be in 
direct relation to the deterioration of the agreement between theory and experiments in 
the determination of the work done by the fluid flow on the Brownian particle. This can 
be understood in a qualitative way. In the presence of a harmonic trapping potential, 
the particle movement along x, y and z are independent. For example, assuming the 
particle to be at x(t = 0) = the protocol acts on the particle trajectory independently 
from the value of z(t = 0). However, if there is a rotational component of the force field, 
a particle will undergo the effect of both the protocol and the rotational force field, 
which, e.g., may result in a longer displacement for z(t = 0) > and in a shorter one 
for z(t = 0) < 0. Clearly, this effect becomes evident only for times long enough for the 
effect of the rotational component to be significant. This explains the widening of the 
experimental work distribution in the case of Ps and Me particles for long times (figure 

ED. 

In principle, the partial differential equation ruling the time evolution of the 
probability distribution function of the work could be obtained also in the presence 
of nonconservative forces following the procedure outlined in [19]. However, while such 
an equation can be explicitly solved in the case of a harmonic potential, leading to 
eq. (131), this is not as easy for the case of nonconservative forces, because it involves a 
nonlinear problem. The discrepancy between the experimental data and the theoretical 
prediction eq. (j3J) observed in figure [3] can thus be attributed to the fact that eq. (j2J), 
and thus eq. ([3]), only takes into account the contribution of the conservative forces to 
the work. 



6. Conclusions and outlook 



We have shown how the presence of rotational radiation forces due to a weak optical trap 
affects the measurement of the work distributions on a driven particle. As we had already 
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observed in [33] , the rotational radiation forces become significative only for a very low 

trapping power, i.e. below 1 mW. This is the regime in which various experiments to 

test thermodynamics relations on colloidal systems have been performed [181 tlH] • We 

have now shown how such spurious effects increase with an increasing refractive index 

difference between the particle and the medium and that they become undetectable in 

the case of a Si particle trapped in water even for an extremely low trapping power. 

These observations will prove useful to minimize the influence of rotational force fields 

in future experiments. 
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x position (nm) 

Figure 1. Schematic of the experiment: (a) a mesoscopic particle performs Brownian 
motion in the harmonic potential well generated by an optical tweezers; (b) at time 
t = 5 s the action of a fluid flow exerts a force on the particle and displaces it from 
its equilibrium position, (c) A single trajectory (grey) and the average over 400 
trajectories (black); the measured average displacement resulted equal to 10 ± 6 nm. 
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Figure 2. Stability of the experimental setup. Experimental power spectral densities, 
which are calculated from the measured particle trajectories |35j for the x (blue squares) 
and z (red diamonds) coordinates, are well fitted by Lorentzians functions (solid lines) 
in the range of interest for our experiments, i.e. from 50 mHz to 500 Hz. 
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Figure 3. Work distributions. Work performed on a silica (first row), polystyrene 
(second row) and melamine (third row) Brownian particle driven by a fluid flow for 
three different times from the beginning of the manipulation protocol. Experimental 
data (histograms) and the theoretical expectation (lines). Deviations are present for 
long times in the case of polystyrene and melamine. 
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Figure 4. Rotational radiation forces. Cross-correlation functions for the p and z 
coordinates of a silica (first row), polystyrene (second row) and melamine (third row) 
optically trapped Brownian particle with the fluid at rest (left) or in the presence of a 
fluid flow (right). 



